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Abstract 

We study the Lagrange formalism of the (rational) Ruijsenaars-Schneider (RS) system, 
both in discrete time as well as in continuous time, as a second example of a Lagrange 1-form 
structure in the sense of the recent paper [23]. The discrete-time model of the RS system 
was established some time ago arising via an Ansatz of a Lax pair, and was shown to lead 
to an exactly integrable correspondence (multivalued map) [15 . In this paper we present the 
full solution based on the commutativity of the discrete-time flows. The closure relation can 
be established by using the equations of motion. The connection with the KP systems is 
also studied. Performing successive continuum limits on the RS system, we establish the 
Lagrange 1-form structure for the corresponding continuum case of the RS model. 

1 Introduction 

The Ruijsenaars-Schneider (RS) system [2D][2T], i.e., the relativistic version of the Calogero- 
Moser (CM) system, is integrable both in the classical and quantum regimes. The classical 
model was discovered in [20] by considering the Poincare Poisson algebra associated with sine- 
Gordon solitons, and was motivated by the discovery in the late 1970s of explicit soliton-type 
S-matrices for some relativistic two-dimensional quantum field theories (such as the massive 
Thirring model, the quantum sine-Gordon theory and the O(N) cr-models). For reasons 
elucidated below, we are interested in Lagrangian aspects of the RS model, which to our 
knowledge have never been considered. An obvious reason for this is that the Hamiltonian 
description corresponding to the system is no longer of Newtonian form, and hence the usual 
connection between the Hamiltonian and the Lagrangian description through the Legendre 
transformation seems no longer suitable. At the same time we are interested in the integrable 
time-discrete version of the RS system, which was proposed and studied in [TE\, where the 
Lagrangian description is more natural than the Hamiltonian one, because the finite-step 
time-iterate can be naturally viewed as a canonical transformation where the Lagrangian 
plays the role of its generating function. In [15) the corresponding discrete-time Lagrangian 
was found, but the continuum limits were not considered so far. As we shall show, the latter 
can be used to derive a natural Lagrangian description for the continuous RS model as well, 
but in the context of what we call a Lagrangian 1-form structure. We will now explain what 
we mean by this latter notion. 
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Recently, a novel point of view was developed on the role of the Lagrangian structure 
in integrable systems, cf. [7J, where it was proposed that the fundamental property of 
multidimensional consistency can be made manifest in the Lagrangians by thinking of the 
latter as components of a difference (or differential) "Lagrange-form" when the flows are 
embedded in a multidimensional space-time. A new variational principle was formulated 
which involves not only variations with respect to the dependent variables of the theory, 
but also with respect to the geometry in the space of independent (discrete or continuous) 
variables. In [7J, this was laid out in the case of two-dimensional lattice equations, whilst in [6 
it was extended also to the case of the 3-dimensional bilinear KP equation (Hirota's equation). 
Furthermore, in [2 3) a universal Lagrangian structure was established for quadrilateral affine- 
linear lattice equations as well as for their corresponding continuous counterparts, the so- 
called generating PDEs of the system. The key property in all these systems, in which in a 
sense the integrability of the system resides, is that the Lagrangian form is closed on solutions 
of the equations of the motion (but not identically closed for arbitrary field values) . This is 
believed to encode the multidimensional consistency of the system under consideration, but 
a formal proof of that connection remains to be given. 

In the case of integrable systems of ODEs, like the equations of motion of integrable 
many-body systems, the Lagrangian form structure is that of Lagrange 1-forms. Recently, 
in collaboration with S Lobb, the authors studied a first example of such a Lagrange 1- 
form structure, namely the case of the discrete-time (rational) Calogero-Moser (CM) sys- 
tem, [231 HI]- The multidimensional consistency of the system in this case is represented 
by the co-existence of two or more independent commuting discrete-time flows in the case 
of three or more particles. Starting with the discrete-time case, we furthermore established 
the Lagrange structure also of the corresponding continuous case by performing systematic 
continuum limits on the discrete-time equations and Lagrangians. Of course, these systems 
exhibit also multi-Hamiltonian structure, in the sense that the classical Hamiltonians in 
involution with respect to a canonical Poisson structure, each of which generates its own 
time-flow in a corresponding time- variable, However, it is not the case that one can perform 
a naive Legendre transformation on each of these Hamiltonians separately to yield a proper 
Lagrangian structure that makes sense as a coherent system. In fact, the higher-order La- 
grangians emerging from such a naive approach would yield rather complicated algebraic 
expressions which seem not to be suitable for further study. However, as we have shown 
in |24j . a proper Lagrangian 1-form structure can be defined for the CM system, in which 
the components of the form are mixed, but polynomial in the time-derivatives, Lagrangians, 
which obey the crucial closure property expressing the commutativity of the flows on solu- 
tions of the equations of the motion. To derive these Lagrangians, the connection between 
the semi-discrete KP equation and the discrete-time CM system, which arises as the pole- 
reduction of the former, was instrumental in order to guide the proper choice of higher-order 
continuum limits obtained by systematic expansions performed on the discrete-time model, 
thus leading to the Lagrangians in the continuum case. Unfortunately, we do not know at 
this stage a Lagrangian of the semi-discrete KP equation in the form used, which would 
perhaps have allowed us to do the pole-reduction on the Lagrangians directly. 

In the present paper, we proceed in the same flavor as with the paper |24j . to establish the 
Lagrange 1-form structure of the discrete-time rational RS system. Whereas in the case of 
the discrete-time CM system there is a direct connection between the Lax matrices and the 
relevant Lagrangians, and where the mentioned closure relation is a direct consequence of the 
zero-curvature condition, such a direct connection is absent in the RS case. Nevertheless, 
a Lagrange structure for discrete-time RS system as was given in [IS], and we show in 
the present paper that this Lagrangian structure (upon a small modification, and when 
extended into a space of multiple discrete times) still exhibits the closure property. Thus, 
this RS case seems to confirm that the Lagrangian 1-form structure is a general feature of 
integrable models in the 1-form case. Furthermore, whereas the discrete-time CM system 
arises as pole-reduction of a semi-discrete KP equation (with one continuous and two discrete 
independent variables) the discrete-time RS is connected by an analogous reduction to the 
fully discrete KP equation (with three discrete independent variables). Thus, the RS case 
seems to be quite a bit richer than the CM case of [23] . It is also interesting to note that one 
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of the lattice parameters of the KP equation (i.e., parameters which are associated with the 
grid size in each direction of the lattice) seems to play the role of the relativistic parameters 
of the RS model. In other words, the non-relativistic limit of the model corresponds to a 
continuum limit in one of the directions of the KP lattice. 

Our focus in this paper is on the rational case of the RS model, even though some of our 
results on the Lagrangians can be extended straightforwardly to the trigonometric/hyperbolic 
case, because, as in the case of the CM system, we prefer all the statements we make to be 
backed up by the explicit solution of the equations of motion that can be constructed in 
this case. For instance, an important ingredient in the structure is what we call "constraint 
equations" , which in addition to the one-dimensional equations of the motion can be verified 
explicitly for the solution obtained. These constraint equations involve the dynamics in 
two discrete variables (corresponding to trajectories in in the space of independent variables 
which involve corners or zig-zags). Since the starting point in the present paper is an Ansatz 
of a Lax pair, rather than a reduction from a KP system (the connection with the lattice KP 
is established a posteriori) the constraint equations emerge from zero-curvature conditions 
on the Lax matrices involving different time-directions. Thus, to establish that the whole 
system of discrete-time equations and constraints has a nontrivial family of solutions backs 
up the role the latter play in the structure. 

The organization of the paper is as follows. In Section 2, we present the full solution 
of the discrete-time RS system, which requires commuting flows in different discrete-time 
directions. In Section 3, the Lagrangian 1-form structure of the discrete time RS system 
will be studied. The closure relation can be derived by direct computation involving the 
equations of motion as well as the constraints. In Section 4, the "skew" continuum limit 
will be derived guided by the exact solution of the discrete-time RS system system yielding 
what we call the semi-continuous RS system. In fact, the latter acts as a generating system 
for the continuum RS system, which allows us in Section 5 to derive the full limit, by 
which we recover the continuous-time RS hierarchy, together with the relevant continuous 
Lagrangians which depend in a mixed way on the various time-derivatives. It is the latter 
system of Lagrangians that constitutes the Lagrangian 1-form structure, and which exhibit 
the closure relation on the solutions of the system. In Section 6, the connection to KP 
systems is presented, by introducing a slight (parameter-)generalization of the original Lax 
representation, characteristic polynomial associated with the exact solution amounting to 
the corresponding lattice KP r-function. Finally in Section 7 we briefly discuss the non- 
relativistic limit, demonstrating that in that limit the relevant lattice parameter effectively 
plays the role of the reciprocal of the speed of light. We end the paper with a discussion of 
some important points on the physical nature of the models studied here. 

2 The discrete-time Ruijsenaars-Schneider system and 
commuting flows 

In this Section we review the discrete-time RS system which has been introduced in |15j . 
This gives us an occasion to introduce appropriate notation which we will use throughout 
the paper. Furthermore, we derive the exact solution of the discrete equations of the motion 
and identify the constraint relations on commuting flows that can coexist in the system. 



2.1 The single-flow RS system 

Following |15) we start with an Ansatz for a Lax pair which in the rational cases takes the 
form: 

= M K 0, L K 0-C0, (2.1a) 



3 



for a vector function </> and an eigenvalue £, in which the matrices L K and M K are given by 

L K = — + L , (2.1b) 

K 

hh T 

M K = + M , (2.1c) 

K 

where 

L °=lt .^^ i . ^d Mo = f] ^±—E i:j . (2.2) 

. . i Xi Xj i ^> i Xi Xj ~r A 

In (|2.ip the are the particle positions, whilst the hi are auxiliary variables which will 
be determined later. The tilde is a shorthand notation for the discrete-time shift, i.e. for 
Xi(n, m) = Xi, and we write Xi(n + l,m) = 2?i, and Xi(n — l,m) = x^, where to is another 

discrete-time variable which will be specified later. The variable k is the additional spectral 
parameter, whereas A is a parameter of the system related to the non-relativistic limit. The 
matrix Eij are the standard elementary matrices whose entries are given by (Eij)ki — SikSji- 

The compatibility condition of the system (|2.1b|) : 



hh T ~ \ hh T \ hh T n/r \f hhT r \ 

I — +L Q \ I — + M ) = I — +M„ f — + L j (2.3) 

From the coefficients of 1/k 2 we derive the conservation law trL K = trL K leading to 

N N 

i=i j=i 
and furthermore, the coefficients of 1/k give 

L hh T + hh T M = M hh T + hh T L , (2.5) 
and the rest produces the equation 

L M = M L . (2.6) 
(|2.5[) and (|2.6[) produce the relations 

N { h 2 h 2 \ N ( h 2 h 2 \ 

Y ^ b. = V -1 - -1 , (2.7) 

y Xi — Xj + X Xi — Xj + X J yxj — xi + X Xj—xi + Xj' 

for all i,j = 1, 2, N. Consequently, both sides of (|2.7p must be independent of the external 
particle label. Thus, we find a coupled system of equations in terms of the variables hi, and 
Xi in the form 



h 2 h 2 \ = 

^— J ' Xi — Xj + A Xi — Xj + X J 



3 = 1 



" ( h) h 2 \ 

\ Xj — x\ -\- X Xj — xi + X J 



where p does not carry a particle label (but would still be a function of " n " and " to "). 
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Lemma 2.0.1. Lagrange interpolation formula: Consider 2N noncoinciding complex 
numbers Xk and y^, where k = 1, 2, N. Then the following formula holds true: 



M « - ») (f - ») n^(» - W) ' 



We now have 



= 1 + 2^ ( x ._ Vi \ t1 n ; :> » = i,...,JV, (2.io) 



which is obtained by substituting £ — xi into (|2.9|) . 

In order to derive the closed set of equations of motion for the variables xi we have to 
determine the variables hi with the use of the Lagrange interpolation formula yielding 



njLiO; + X)(xi - xj - A) 
-P " N , STTjV , . ( 2 - lla ) 



- acj) Uj=i( x i - Xj) 



r 2 0&i - x,j + X)(xi — Xj — A) 

h i = p-^t^z — „ (2.11b) 

llj=£i\ X i X j ) LLj=l\ x i X j) 

for i — 1, 2, N which we obtain the following system of equations 

P tt (xj — Xj + A) _ -pj- (xj — x J )(x l — Xj + A) 
P 1 = 1 fa - - A) 1 = 1 fa - x,-)fa - 2,- - A) ' 

(|2.12[) can be considered to be the product version of (|2.1ip . is a system of N equations for 
N + 1 unknowns, x\, ...,xn and p. There is no equation for p separately, and thus it should 
be a priori given in order to get a closed set of equations. If p is a constant implying pj p 

to be equal to unity, we obtain the equations of motion of what we would like to call the 
discrete-time RS system corresponding to the " ~ " direction. 

The exact solution of the equations of motion can be derived in a way similar to the 
continuous case of the rational RS model, cf. .15: , cf. also [20l [19] using the explicit form of 
the Lax matrices (|2.1b[) . The details are given in Appendix A, and we can summarize the 
conclusion as follows. 

Proposition 2.1. Let the N x N matrix function of the discrete variable n, Y(n), be given 
by 

Y(n) = (pi + A)- (V(0) - (pi + A)" , (2.13) 

subject to the following condition on the initial value matrix 

[y(0),A] = AA + rank 1 . (2.14) 

then eigenvalues Xi(n) of the matrix Y(n) coincide with the solutions for particle position of 
the discretetime RS system, i.e. they solve the discrete equations of motion (|2.12D . 

As a corollary, we have that equivalently the solutions can be found from the secular 
problem of the matrix: 

Y(0)-npX(pI + L (0))~ 1 . (2.15) 
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where 



2.2 Commuting discrete flows 

Following the construction in [24] we introduce another rational Lax pair in the form 

4> = N K 4>, K K 4> = C4>, (2.16a) 
for the same vector function <f>, and 

K K = — + K , (2.16b) 

K 

N K = — + JV , (2.16c) 

N h h N h h 

= E 7r^i E " • " ld N ° - E ■ (2 - 16d) 

z 5i 7 — 1 J ij — 1 

which provides the flow in an additional discrete-time variable m. Similar to the Lax pair 
in (|2.1j) . in (|7.2p the ki are auxiliary variables which will be determined later. The hat 
is a shorthand notation for the discrete-time shift, i.e. for Xj(n, m) — Xi, and we write 
Xi(n, m + 1) = Xi, and Xj(n, m — 1) = Xj. 

Obviously, the compatibility relations for (|2.16[) can be analysed in a very similar man- 
ner as to the ones for (|2.ip . Thus, we find a coupled system of equations in terms of the 
variables fcj, and Xi in the form 

^ ^ ^ Xi — Xj "I - A Xi — Xj -\- A y 



where g does not carry a particle label. The system (|2.17[) can be resolved once again by 
using the Lagrange interpolation formula (12.101) yielding the resolution: 

, 2 Yl'j ! =1 (x l - Xj + \)(xi ~ Xj - X) /„ no \ 

K = _ 9^Hv jv ~ ' (2.18a) 

n^o^ _ x j) Y[j=ii x i ~ x j) 

% = g n ^ i(£ r_ a: i +A) f < :^" A) , 

rij^c^i ~ x j)Y[j=i( x i ~ x j) 

for i = 1, 2, JV, and from which we obtain the following system of equations 

q -pr (Xj - Xj + A) = -pr (X. t -Xj){Xj -Xj + A) _^ 
9 fj[ ( x i - x j ~ A ) } = i 0* - - % - A ) ' 

The product version (|2.19l) of (|2.18p . thus yields a system of JV equations for JV + 1 unknowns, 
xi, ...,xjv and g. There is again no equation for q separately, and thus it should be a priori 
given in order to get a closed set of equations. If q is a constant implying q/q to be equal 
to unity, we obtain the equations of motion of the discrete-time RS system corresponding to 
the " ~ " direction in terms of the discrete-time variable m. Thus far, so similar. 

Assuming now that the dependent variables depend simultaneously on both discrete time 
variables n and m, then to obtain a univalent solution of the equations of motion we must 
require that both flows, in the " ~ " direction and the " ~ " direction, commute. If so, then 
we can fix a value for n and solve the equations in the " ~ " direction similarly as before, 
leading to the matrix Y which depends on n and m as follows: 

Y(n, m) = (ql + A)-™ (V(n, 0) - ^f^) {ql + A) m , (2.20) 
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subject to the constraint on the matrix Y at m = 0: 



[Y(n, 0), A] = AA + rank 1 . (2.21) 

In order for this scenario to work there must be further constraints on the flows. This 
will lead to a system of "constraints" which can be readily obtained from the compatibility 
between Lax pairs (|2.1c[) and (|2. 16c|) . In fact, the compatibility relations between the simul- 
taneous eigenvalue problems for the matrices L K and K K produces the following relations 
(see appendix iBj) 

h\ = (3k 2 , (2.22) 



where 

2^.7 = 1 K i 

Using (l2~TTj) and ([2~T8)l . we obtain 



T h 2 

I 

N 



— = f] — ~ % - X ^ Xi ~ ^ , (2 24a) 

q/3 jj^ (xi - xj - X)(xi - Xj) ' 

P_ _ TT - Xj + X)(xj - Xj) 

q/3 ( 2 ; i _ + A) (a;, - Xj) ' 

We will refer to relations (|2.24al) and (|2.24bp as the constraint equations which guarantee the 
commutativity between the discrete-time flows with shifts " ~ " and " ^ " in the variables n 
and m respectively. 



Proposition: The eigenvalues X\{n, m), . . . , xm{ti, m) of the N x N matrix 

Y(n,m) = { P I + A)- n (qI + A)- m Y(0,0){ P I + A) n (qI + A) m 

-npXipl + A)" 1 - mqX{qI + A)" 1 (2.25a) 

in which the initial value matrix Y(0, 0) is subject to the condition 

[Y(0, 0) , A] = AA + rank 1 , (2.25b) 

obey both the discrete-time Ruijsenaars- Schneider systems given by eqs. (I2.12j) and (|2.19[) as 
well as the systems of constraint equations given by (|2.24aj) and (|2.24bp . 

In order to make a connection with an initial value problem, we mention that the initial 
value matrix Y (0, 0) can be obtained from the diagonal matrix of initial values X(0, 0) by a 
similarity transformation with a matrix 17(0,0) which is an invertible matrix diagonalizing 
the initial Lax matrices Z(0,0) and K(0, 0). To find the latter, we need the initial values 
Xi(0, 0), Xi(l, 0) and 2^(0, 1), {i = 1, . . . , N). We note that the secular problem can, hence, 
be reformulated as one for the following matrix 

X(n,m) = X(0,0) -npAL _1 (0,0) - mqXK^iO, 0) , (2.26) 

and hence the solution is provided by the roots of the characteristic equation: 

N 

Px(x) = det(xl - X(n,m)) = J[(x - x l (n,m)) . (2.27) 

i=i 

Remark: We would like to mention at the end that if p and q are no longer constants the 
compatibility between (|2.13[) and (|2.20l) produces the conditions 

pq = qp and p + q = q + p . 

From these two equations, the only possible answers would be p — p{n) and q = q(m) 
implying that p and q must be functions of their own a discrete variable. 
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3 The Lagrangian 1-form and the closure relation 



In the CM case [24] , we obtained Lagrangians 1-form structure through the connection of the 
Lax representation. Here we also have the Lagrangian 1-form structure for the RS system, 
but the establishment is more difficult as connection through the Lax representation is no 
longer relevant. In this Section, we will first derive the Lagrangian 1-form for the discrete- 
time RS system and then establish the closure relation. We begin with rewriting (I2.12j) and 
(j2~T9|) in the forms 



JY 



(hi(xi — Xj + A) — m(x, — Xj — A)) = 



3=1 



JY 



N 

£(ln(a 

3=1 



,•) + ln(xj - Xj + A) 



ln(xj - Xj) - \tl(x 1 - Xj - A)) , (3.1) 



(ln(xi — Xj + A) — \n(xi — Xj — A)) — 



3 = 1 

3±i 



(ln(xi — Xj) + hi(xi — Xj + A) 

3 = 1 



— ln(xi — Xj) — ln(xi — Xj — X)) . (3.2) 
It is easy to show that the actions corresponding to these equations of motion are given by 

S (n) = ^2^(n)(^(n),x(n + l)) , 



(3.3) 
(3.4) 



where 



N 1 N 

"^(n) = tt( Xl ~ ~ ~ %3 ~ A )) - n (f( X i ~ X 3 + A ) 



(m) 



with 5 = J] 



AT 



i,3=l 



i,3 = l 



- + A)) - In 



N N 
Y (f( X i ~ %) _ f( X i ~ X J ~ A )) ~ o ~ X 3 + A ) 



i,3=l 

+ - Xj + A)) - In 



i,3 = l 



(3.5) 



(3.6) 



the function f(x) given by /(x) = xln(x). The discrete-time Euler- 



Lagrange equations read 



as?, 



(n) 



(™) 



<9x,; 



<9x,; 



, and 



(m) 



(m) 



<9x,; 



<9x,; 







which lead to (I3.1[) and (I3.2j) . respectively. 

The additional terms in (|3. 51) and p. 61) . containing the differences of the centre of mass, 
are needed in order to account for the constraint equations (|2.24aj) and (|2.24b|) as they arrive 
from the EL equations on discrete curves, which is a connected collection of line segments 
(i.e. elementary links on the lattice) with or without end points (i.e. closed or non-closed), 
involving corners (vertices connecting line segments with different directions). 



Theorem 3.1. (The closure relation) 

_Sf(„) (as, i) - Jz?( n ) (as, x) - _Sf (m ) (as, x) + «5f (m) (as, x) = , 



(3.7) 



holds on the equations of motion and the constraint equations. 
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Proof. (|3.7[) can be written in the form 













A 




Xj 




Xi - 


Xj 




Xj — 


A 


Xi 


— Xj 






— Xj 




Xi 


Xj — 


A 


Xi - 


Xj — 


A 


X i 


— Xj 




Xi 


Xj — 


A 


Xi 


— Xj 




Xi 


— Xj 




Xi 


Xj — 


A 



In 



In 



Xi — Xj -\- A 



X i X j A 



X> i X ^ I 



X> £ Xj j 

Xj In 



X i X j A X j 



cc) - if (n) (x, cc) - 5?( m )(x, x) + 3?( m )(x, cc) 

N ~ ( 

= Xj I In 

i,J=l V 

w ~ / 

i,J=l V 

N (~ 

lnjq^ - 

N / 

Using (EH]), (EH]), (I2.24aj) and (|2.24b|) . we have 

_Sf (n) (cc, cc) - _Sf (n) (cc, cc) - _Sf (m) (cc, cc) + _Sf (m) (cc, cc) 
= y^(ccj + ccj - Xi - Xi) In 



X j ^"i *^ J 



In 



Xi 


Xj 


- A Xi 


Xj 


- A 


-In 


Xi 


Xj 


\- A 


) 




— Xj 


-\Xi 


— Xj 


- A 


Xi 


— Xj - 


1- A 





p_ 

q0 



ln\qy/ji 



N 



-In 


p 
































Xi 


Xj — 


A a;, - 


Xj 


- A 


-In 


Xi 


Xj 


1- A 


) 


Xi 


Xj — 


A x i ~ 


Xj 


-A 


Xi 


— Xj - 


1- A 





(3.8) 



(3.9) 



Using the fact that the last line of (|3.9[) vanishes on the exact solution (|2.25a[) and 5 — 5 
H + H = 0, then we have 



^(n)(x,x) -J? (n) (cc,cc) -«Sf (m) (cc,cc) +^f (m) (cc,cc) = 



(3.10) 

□ 



n ; 



it i 



V 



71: 



Figure 1: Deformation of the discrete curve T. 

In [53] we described what we mean by the Lagrangian 1-form, but let us reiterate this 
here for the sake of self-contained of this paper. Let represent the unit vector in the lattice 
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direction labeled by i and let any position in the lattice be identified by the vector n, so 
that an elementary shift in the lattice can be created by the operation n i— > n + e^. Since 
the Lagrangian depends on x and its elementary shift in one discrete direction, it can be 
associated with an oriented vector ej on a curve Ti(n) = (n, n + ej), and we can treat these 
Lagrangians as defining a discrete 1-form Jz?j(n) 

ifi(n) = S?i(x(n),x(n + ei)), (3.11) 

which satisfies the following relation 

J£i(x{n + ej), x(n + e, + e^)) — Jz?j(a;(n), x(n + ej)) 

-^fj(x(n + e t ),x(n + ej + ej)) + 3?j(x(n),x(n + ej)) = . (3.12) 

Equation p,12[) represents the closure relation of the Lagrangian 1-form for the RS system 
and it can be explicitly shown holding on the level of the equations of motion, and as well 
as constraints. 

Choosing a discrete curve T consisting of connected elements we can define an action 
on the curve by summing up the contributions Jzfj from each of the oriented links I\ in the 
curve, to get 

S{x(n);T) = ]T ^: (x(n), x(n + e^)). (3.13) 

The closure relation (|3.12l) is actually equivalent to the invariance of the action under local 
deformations of the curve. To see this, suppose we have an action S evaluated on a curve 
F, and we deform this (keeping end points fixed) to get a curve V on which an action S' is 
evaluated, such as in Figured! 

Then S' is related to S by the following: 

S' = S — J^i(x(n + ej),x(n + e, + ej)) + JS?i(x(n),x(n + e^)) 

+J£j(x(n + ei),x(n + ej + e.;)) — Jg?j(x(n), x(n + ej)). (3-14) 

Equation (|3.14[) shows that the independence of the action under such a deformation is locally 
equivalent to the closure relation. The invariance of the action under the local deformation 
is a crucial aspect of the underlying variational principle. In Appendix [CJ we demonstrate 
how to derive the discrete Euler-Lagrange equation the variational for some specific discrete 
curves. 

4 The semi-continuous limit: The skew limit 

In this Section, we study a continuum analogue of a previous construction in Section [5] by 
considering a particular semi-continuous limit. Since the exact solution (I2.25ap contains two 
discrete variables n and m, we could perform a continuum limit on one of these variables 
separately, while leaving the other discrete variable intact, and thus obtain a semi-continuous 
equation with one remaining discrete and two continuous independent variables. Alterna- 
tively, we can first perform a change of independent variables on the lattice and subsequently 
perform the limit on one of the new variables. The advantage of the latter approach over 
the former is that it often leads in a more direct way to a hierarchy of higher order flows. 
Adopting the latter approach in this section, we use a new discrete variable N := n + m, and 
perform the transformation on the dependent variables by setting x(n, m) <-> x(N,m) =: x, 
which leads to the following expressions for the shifted variables: 

x = x(n + l,m) i y x(N + l,m)=:x, 
x = x(n, m + 1) n- x(N + 1, m + 1) =: x , 

x = x(n + l,m + 1) i — ^ x(N + 2, m + 1) =: x . 



10 



Rearranging the terms in (|2.25a[) . we have 

' ql + A 



Y(N,m) = { P I + A)- N 

P 



pi + A 



Y(0,0) 



pI + A 



pI + A ql + A 



(pi + A) 



(ql + A 
\pl + A 



(4.1) 



We perform the limit n — > — oo, m — > oo, e — >■ while keeping N fixed and setting e = p — q, 
such that em = r remains finite. Focusing on the penultimate factor in (14.11) we have that 



lim 1 - 



pI + A / 

so that the exact solution takes the form 



= lim 1 — 



m{pl + A] 



Y(N,t) = ( P I + A)- N eTT 



Y(0,0)- 



r\A 



NpX 

pi + A + (pi + A) 2 



= e p'+ A , 



(pl + A) N e~^4 



(4.2) 



(4.3) 



This equation represents the full solution after taking the skew limit. The position of the 
particles Xj(N,r) can be determined by computing the eigenvalues of (|4.3[) . 

4.1 The skew limit on equations of motion and constraints 

We rewrite the equations of motion (|3.2j) in terms of the variables (N,m) 

N N 

^ Qn(xj - Xj + A) - ln(xj - Xj - A)) = ^ f ln(x 4 - xj) - ln(x. ( - xj + A) 



3=1 



3=1 



ln(xi - x,-) - ln(x 4 — Xj — A) J , (4.4) 



Using a Taylor expansion 



^ , . . . dx e 2 d 2 x 

x = x(r + e)=x{r)+e- + -— 2 

. . . . dx e 2 d 2 x 



(4.5) 
(4.6) 



and collecting terms in order 0(e°), we have the equations of motion for the RS system 
corresponding to the " N " variable 



N 



N 



^ (ln(xi - xj + A) - ln(xj - xj - A)) = (Info - xj) - Info - Xj + A) 

3=1 



3 = 1 



+ ln fo - 5j) - ln fo ^ x j - A )) 



and 0(e), we have 



E 

3=1 



1 



dx 3 

<9r \ Xi — x 



A 



9x„- 



9t y Xj — Xj + A Xi — Xj 



= 



(4.7) 



(4.8) 



which are the equations of motion for the RS system corresponding to the r variable. 

Similarly, changing the variables to (N, r) in (|2.24a[) and (I2.24b[) and collecting terms in 
order 0(e), we have 



N a~ 

Eoxj 



^-^ dr V x.i — xa — A Xi — Xj 

j=i J 

1 ^ 9xj / 1 1 



3=1 



c?r V Xj — x,- + A Xi — x, 



(4.9) 
(4.10) 



11 



The summation of these two yields (|4.8[) . while the difference gives 



N 

E 



x, — A 



9t \ Xi — Xj + A x j — Xj 



= 0, 



(4.11) 



which are the constraint equations after taking the skew limit. 
4.2 The skew limit on action 

In this Section, we proceed exactly the same steps what we did in [2~4"| . First, we observe 
that eq. (|4.7[) can be once again be obtained by implementing the usual variational principle 
on the following action S(n) given by 



S, 



(N) 



N N 

E = E E (/( x * - **) - - % - A )) - o E /( x * - x * + A ) 



N \ i ,j=l 



9 E ^ Xi ~ x j + A ) - m ^ E( Xi ~ Xi ) 



(4.12) 



i = l 



where now the Lagrangian -S?(n) involves variables shifted in the discrete variable N instead 
of the original variable n, and the corresponding discrete Euler-Lagrange equation reads: 



dxi \ dxi 



= 0, 



(4.13) 



yielding (|4~71) . 

Second, we observe that eq. (|4.8p can be once again be obtained by implementing the 
usual variational principle on the following action S7 T ) given by 



S(r) = J J dT^ (T) (x(N -l,r), 



ax(N ,r; 
8t 



(4.14) 



which is obtained by taking the skew limit together with anti- Taylor expansion of (|3.4[) and 



JY 



Ok 



Jz? (t) = ( ~jfp ( ln l Xl ~ x .j ~ A l ~ ln l x * _ x j 



1 N 



2 ^ V dr 

iV 



^ (In ft - ^ + A| - m ft - Xj - - A|) + ^ - ^ 



E u( x *- x ») 



<9r 



ln 



(4.15) 



The Euler Lagrange equations 



d&( T ) _ d_ f OJ^r) 

9xj dr \d(dxi/dr 



yield KEl . 



5 The full limit 



= o 



(4.16) 



In the previous Section, we took the continuum limit on the discrete variable m, leading 
to a system of differential-difference equations. The full continuum limit, performed on the 
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remaining discrete variable N, will lead to a coupled system of poles in the first instance, from 
which a hierarchy of ODEs can be retrieved, which is the RS hierarchy. How to perform 
this limit is inspired by the structure of the solutions of (I4.3[) . Performing the following 
computation, 

Y(N,r) = (1+^ eM I+ Try(0,0)e-f( I+ Tr (i + j " 

( A\ _1 t\A / A 
-NA 1 + - +—5- 1+- 

V PJ P V P 

= e -Nln(l+^) + j(l+f)- 1 Y(0j0)e Nln(l + A)-r( 1+ A)^ 

V pJ p V P/ 



We now introduce 

r N 2r N 3t N 

tl — H > 12 — r 5" ' *3 — — 7 H — 7 7 ; (5.2a) 

p pa pz p+ p.i 

and expand (|5.1[) with respect to variable p. We have 

Y(ti,t 2 ,t 3 ,...,N) = e- Atl+A2 ^- A3 *+-Y(0,0,...)e Atl - A2i ?+ A3 *+- 

-NA + AAti + A 2 At 2 + A 3 At 3 + ... , (5.3) 

which is a function of time variables (tx, N). The positions of the particles Xj (ti, ti, t$, 

can be computed by looking for the eigenvalues of (|5.3|) . The explicit expression of the so- 
lution for the RS can be obtained from the secular problem for the matrix 

X(0, 0) - £ + £(0, 0)Ati + i 2 (0, 0)\t 2 + L 3 (0, 0)Xt 3 , (5.4) 

where £ = NA and X(0, 0) = *7 _1 (0, 0)F(0, 0)17(0, 0) and L(0, 0) = E/ _1 (0, 0)AJ7(0, 0). The 
solution (|5.3I) involves TV-time flows for the RS system. The next solutions in the hierarchy 
can be generated by pushing further on with the expansion. 

5.1 The full limit on the equations of motion 

We now would like to see what would result from taking the limit on the equations of motion 
First, we introduce 

dxi dXi dh , dXt dt 2 , dX, dt 3 



and 



dr dti dr dt 2 dr dt 3 dr 

1 dX, 2 dX t 3 dX t 
H + ... 

p 2 dt\ p 3 dt 2 p A dt 3 



, IdX, 1 (\d 2 X l dXA 1 fldXi d 2 X { dX t 
Xj - A + 



and 



x.i + A 



(5.5) 



p dh p 2 \2 dt\ dt 2 J p 3 \6 dt 3 dhdt 2 dt 3 

i/ia 4 x, i d 3 x t \d 2 x l d 2 x t \ 5 , 

? U^r " 2^ + 2^r + + 0{1/p } ' (5 - 6) 

ldXi 1 /la 2 ^ <9JQ\ 1 / ldXi d 2 X l dXi 



p dti p 2 \2 dt\ dt 2 j p 3 V 6 dt 3 dt x dt 2 dt 3 

V U^f + 2^ + 2^f + > + 0{l/p ] ■ (5J) 
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Then we expand (|4. 8|) with respect to the variable p together with (|5.2a[) . We find that 
♦The leading term of order 0(l/p 3 ) gives us 

1 1 2 



A' 



d 2 X,, I dXj y\ dX 3 



dt\ I dti 



dt\ V X^ — Xj -\- A Xi — Xa — A Xi — Xj 



= 



(5.8) 



which is the equations of the motion for the continuous RS system. 
♦The term of order 0(l/p 4 ) gives us 



<9 2 X- IdXi d 2 XidXi 



dtitil dh dt\ dt 2 

N 



dXj 
dt 2 



l d 2 Xj 
A dt{ 



3=1 

ld 2 X 



dX, 



1 



Xi — Xj -\- A X% — Xj — A X i — Xj 



i 



2 dt\ \Xi - Xj ■. - A X, - X 3 + A 

2 



,1 ^ 

2 V 9t 2 



1 



(X, - X, - A) 2 (X - X, + A)s 



Q 



(5.9) 



This equation represents the next equation of motion beyond the usual continuous RS in 
the hierarchy. We will stop at this equation, but we can actually get the higher terms of 
the equation in which the variable h and higher order time- flows must be taken into account. 



The full limit of (|4.11l) in the order 0(l/p 2 ) gives 

; _ 9 OAi / <).\, 

A dh 



1 



1 



2dXi_dX i /dXi ^dX i 

dt 2 ' dh f-f dh yXi-Xj + X Xi-Xi-X 



(5.10) 



which is the constraint equations for the full limit. 
Using (|5.10[) . we can simplify (|5.9I) into 



d 2 X, OX. 



dh 



JY 



TJ ~dh ^ 

3 = 1 



dXj 
~dh 

1 



1 



Xi — Xj -\- A Xi — Xj — A Xi — Xj 



1 dX, dX, 

2~dh~d~h y(Xi - Xj - A) 2 (Xi - X 3 + A) 2 



= 



(5.11) 



Note that (|5.11[) can be obtained directly from the full limit in order 0(l/p 3 ) from the 
combination of the relations 



1 9xj 

~ p = h 

1 ^ 

n ~ 2^ ~fW 



dr \ x. 



1 



A Xi 



j'=i 



dr \ Xj — xh — X x. 



(5.12) 
(5.13) 



which are the backward shift and forward shift of (|4.9p and (I4.10j) . respectively. 



5.2 The full limit on the action 

We will follow the steps in [24] in order to obtain the full limit on the action. We now take 
the action to be of the form 

5[aj(N,r);r] = P dr-JSf (T) (x(N,r),i(N,r)) +5^Jgf (N) (aB(N,T),a;(N + l,r)) , (5.14) 
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where the first term belongs to the vertical part and the second term belongs to the horizontal 
part of the curve F. 

Using anti- Taylor expansion, the action now becomes 



S[sc(N,r);r] 



drJSf (T )(a;(N,r),x(N,r))+ / dNjSf (N) (aj(N,T),a;(N + l,r)) , (5.15) 



where we do not need to take into account the boundary terms coming from the expansion, 
because they are constant at the end points and do not contribute to the variational process. 
We now perform a change of variables (r, N) i-> (ti 7 t 2 ) by using (|5.2ap 



dr = —p^dt 2 — p 2 dt\ , 
dU = p 2 dt 2 + 2pdh , 

and also expand the Lagrangians with respect to variable p. We obtain 

dX(h,t 2 ) 



(5.16a) 
(5.16b) 



S[X(t u t 2 );T] 



/ dh2 {tl) [X(h,t 2 ), 

[ t2(2 \it * (x(t t \ ^ihiM dx ^) 



where 2^ tl ) and Jz?(t 2 ) are given by 



(*i) 



E^ ln 



dh 



N 



dX 



Y d ^{\n\X i -X j -X\-\n\X i -X j \) 



of which the Eulcr-Lagrangc equation 

fljsf (tl) o ( aj2? (tl) 



OX. 



gives exactly eq. (|5.8[) and 



2, 



(*2) 



N 

E 



dXj 
dt 2 



In 



dXt 



dh 



1 



(5.17) 



(5.18) 



(5.19) 



2A V dh 



dXj 
' dh 



N 

E 

1 dXi dXi 



ay 

--i (In |JQ - ^ - A| - In \Xi - X, 

Ol2 



1 



2 <9ii 9*! Xj - Xj + X 



(5.20) 



We see that the Lagrangian 2u 2 \ contains derivatives with respect to two time flows h and 
t 2 . We observe that the equations of motion (|5.11l) require the Euler-Lagrange equation in 
the form 



82, 



(*2) 



d ( as?. 



(*2) 



^ 9t 2 U(^) 



Furthermore, we find that 
dJ? (t2) 



d m) xdh 



N 



E 



dX j 



= 



1 



1 



2 dXj dXj dx 

dt 2 I dh ' ^ dh V Xi - X* + X Xi- Xi - X 
j=i ■> ■> 



(5.21) 



= , (5.22) 



which is identical to term with the full limit in order 0(l/p 2 ) of the constraints (|5.10l) . 

Here we obtained the hierarchy of Lagrangians for the RS-model through the full contin- 
uum limit. Obviously, higher Lagrangians in the family can be generated by pushing further 
on with the expansion. 
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5.3 The full limit on the closure relation 

We take the full limit on the discrete closure relation 13.71 leading to 

Theorem 5.1. We find that the continuous version of the closure relation between t\ and 

*2 



0JS?r 



(*i 



dh dt 2 

which holds on the equations of motion and constraint 
Proof. : We find that 



and 



9i 2 



N 

E 



a 2 x 



N 

E 



a 2 x. 



In 



ax 



dt 2 dh 



+ 



ax. ax 

'~dh~dh 
dh dt 2 



dti 
[In \X. 
1 



a 2 x 



ati9t 2 
x 



, A| - In |X 
1 



i 



A i — Xj 
1 



X 



X-A 



Xj — Xj 



(5.23) 



(5.24) 



(*2) 



We find that 



<■>//; 



iV 

E 



a 2 x 

a 2 X 



In 



ax 



ax a 2 x / ax 1 ax a 2 x 



+ 



E 

ax aXj 

~dh~dh 
dXj dXj 
dh dh 
1 d 2 Xj dXj 
'2~dtf~~dh 



dhi-^ 



dh dt\ I dh A dh dt\ 



X,--A|-hi|X-Xj] 



X — X; — A X 



X 



X-X 



, - A X 
1 



X 



A^ — Xj -\- X Xj — Xj — A 



i)l. 2 



<*l) _ 



E 



/ax 


A 2 ax 


1 


I dti 


J a<! 


. (X — Xj 




gives 






ax. 


ax 




1 


H - 

ati 


ai 2 


X — Xj — X 



X — X 



JV 

E 



, a 2 x 
at x at s 



(5.25) 



i ax a 2 x a 2 x 



A ati dt{ dhdh 



i a 2 Xj- ax 

'2~dtf~dh 



N 

E 



ax dXj 

' dh dt 2 
1 



ax a 2 x / ax 
at 2 ai 2 / dh 

i 



Xj — Xj 



X{ — Xj -\- A 



Xj — - Xj — A 



(dXf 


\ 2 ax " 


\dh y 


/ dh 



(X - X, - A) 2 (X - Xj + A) 2 



(5.26) 
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Dividing (j5T26l) by P± we find that 



d^ {tl) d^ {t2) = dx t I d 2 x, / dXj d 2 x t dx t I ( dJQV 1 d 2 x, 

dt 2 dti ^ dh y dhh 1 dti dt\ dt 2 ' \ dt 2 J X dt\ 

N 



E 



dX n 



1 



1 



$£2 V X{ — Xj -\- X Xi — Xj — A Xi — Xj 



1 d 2 X n 



2 dt^ \ X h — Xj -\- X X{ — Xj — A 



i fdXi 



i 



i 



2 v dt 2 ) v {Xi - x, - a) 2 {x, - Xj + xy 



(5.27) 



Using (fOT)) becomes 



2^ m, M 



V ^ / d 2 X, dXi I fdX^ 1 



dt 2 dh ^ dh \ dt 2 dt 2 I V dt 2 

N 



-2 > J - 



3=1 

Inserting (|5.8[) . we have now 



1 



1 



(9^2 V Xi — Xj -\~ X Xi — Xj — A Xi — Xj 



(5.28) 



dt 2 dh 



',3 = 



a.v, ( -).v, a.v, ax,- 



^ V dt i dt 2 dh dt 2 J\X i -X j +X Xi-Xj-X Xi-Xj 



The first term of f|5 . 29[) is the antisymmetric function, hence vanishes. We now have 



dt 2 dh 



= 



(5.29) 

(5.30) 
□ 



6 The connection to the lattice KP systems 

In contrast to the CM case [23], where we started with a semi-discrete KP equation, and 
applied a pole-reduction to it to a yield a compatible CM system, here we start from RS 
system and reconnect it to the fully discrete lattice KP systems. In [15], the connection 
between the RS system and the KP system was established for the trigonometric case, but 
here we will focus on the (simpler) rational case as it clarifies the situation more clearly, 

We will develop now a scheme along the lines of the papers [55] [HI HI US] . Starting from 
the "solution matrix" Y(n,m) of (|2.25a[) . we will introduce the relevant r-function as its 
characteristic polynomial: 

r(0 = det(£J - Y) , (6.1) 
where Y = Y(n, m, h) is a function of three discrete variables 

Y = ( p I + A)-^Y(pI + A)--^- , (6.2a) 

pi + A 

Y = (ql + Ar'Yiql + A)--^— , (6.2b) 

ql + A 

Y = {rl + A)~ 1 Y(rI + A) ^— , (6.2c) 

rl + A 
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with r is a lattice parameter corresponding to the " ~ " and we show that r(£) plays the role 
of the tau-function of a discrete soliton system. 

To derive the equations directly from the resolvent of the matrix Y , we proceed as follows. 
First, we perform the simple computation 

f(£) = det(£ + X-Y -rs T ) , 

= dct((e + A-r)(l-? S T (e + A-Y)- 1 )) , 
= T^ + X)(l-s T ^ + X-Y)- 1 r) , 



then wc have 



=l-s T (Z + X-Y)- 1 (p + A)- 1 r = v p (Z + X) . (6.3) 



r(S, + X) 

The reverse fraction of Eq. (|6.3[) can be computed by processing the same computation 

r(0 - det(£ - A - Y + rs T ) , 

= dct((e-A-r)(l+? S T (j-A-F)- 1 )) , 
= f(e-A)(l+5 T (e-A-^)- 1 ?) , 



then wc have 



n ° - l + s T (p + A)- 1 (e-A-F)- 1 ? = w p (e- A) . (6.4) 



f(£-A) 

From (|6.3[) and (|6.4[) . we have the relation 



r(i ' vvy- A: — L- . (6.5) 



r(C-A) ' v p (0 

The same type of the relation for the other discrete direction can be obtained in the forms 

rlii -w,(f-M) = ^, ( 6 - 6a ) 



- w r (e - 77) = -4tt ■ ( 6 - 6b ) 



In order to derive discrete KP equations for r(£), w and v, we introduce the N-component 
vectors 

u«(0 = K-VJ-^o + AJ-V, (6.7a) 
<ti 6 (0 = a^ + A)- 1 ^-^)- 1 , (6.7b) 

as well as the scalar variables 

■MO = s7 > + A r'(? - ^T 1 ^ + A )~ lr • ( 6 - 8 ) 

We now consider ()6.7a[) which can be written in the form 

y-Y)u a (0 = (a + A)- 1 (p + A)?, 
(£- \-Y + rs T )u a {0 = (p-a)(a + A)- 1 r + r , 

K-A-?K(6 = (p-a)(a + A)- 1 ? + ?(l- S T t ta (e)) , 

«a(0 = (p-o)S a (C- A)+v a (OM (e- A) , (6.9) 

with «o(0 = (£ - 

The same process can be applied to (|6.7h>[) and we obtain 

= (P - &) V* + A) + w 6 (0 V)« + A) , (6.10) 
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with *u (£) = s T (£-Y)- 1 . 

Another type of relation can be obtained by multiply s T (b + A) -1 on the left hand side 
of (HH). We have 

^(b + A)- 1 ^^) = (p-a)* r (& + A)- 1 fi (£-A) 
+v Q (05 T (6 + A)- 1 M (e-A) , 
s T ( P + A)(b + A)- 1 u a (C) = (p-a)So6(£-A) + v„(Ow 6 (£-A), 

VaK)w6«-A) = 1 + (P-6)S„6(0-(P-O)5ai($-A). (6.11) 

Similarly, multiplying the right hand side of (|6.10l) . we obtain 

w b (£)v a (£ + A) = 1 + (p - 6)S ab (£ + A) - (p - a)S ab (0 . (6.12) 

By proceeding the similar steps, we can derive the relations in other discrete-time directions, 
namely 

v„(Ow 6 (£-m) = 1 + (« - b)S ab (£) - (q - a)S ab (£ - p) , (6.13a) 
Va(OW&(£-77) = l + (r-6)S a6 (0-(r-a)Sa&(£-77) , (6.13b) 

Using the identity 

w b (£ - A - 7?)v a (C - y) = w fc (g- A-7?)v a (g- A) w b (£- A-p)v a (£-p) 
w b (£ - p - 7?)v a (^ - v) w b{£, - A - p)v a (£ - A) w b (£ - p - 7?)v a (£ - p) 

we can derive 

1 + (p - b)S ab {£ -r))-(p- a)f ab (£ - A - t?) 

1 + (g - b)S ab (£ - j]) - (g - a)S ab (£ - p - 77) 

= 1 + (r - b)^ ab (e - A) - (r - a)f n6 (g — A — 77) 

1 + (g - 6)S ab (£ - A) - (g - a)f ab (£ - A - p) 

1 + (p - b)S ab (£ - p) - (p - a)S ab (£, — A — p) 
x ~ ZE ) 

1 + (r - b)S ab (£ - p) - (r - a)S ab (£ - p - 77) 



(6.14) 



(6.15) 



which is a three-dimensional lattice equation which appeared first (in a slightly different 
form) in [TT]. Effectively, this is the Schwarzian lattice KP equation which in its canonical 
form was first given in [3] , cf . also [35] . 

We now multiply s T on the left hand side of (|6.9I) leading to 

s T «a(£) = (p - a)s T u a (£ - A) + v a (£)s T tt (£ - A) , 
s T (p + AK(0 = (p-a)(l-v a (£-A))+v a (£)S T w (£-A) . (6.16) 

Introducing 

woo(e)=^ T (^^)" 1 r, (6.17) 

(|6.16[) can be written in the form 

(p + 2oo(£ - A))v tt (0 - (p - o)v„(0 = a + s T Au a (£) . (6.18) 

Another two relations related to the " ^ " and " ~ " directions can be automatically obtained 

(g + u 00 (£, - p))v a (£) - (g - a)v a (C - p) = a + /Au„((), (6.19a) 
(r + u O0 (£-77))v a (O-(r-a)v a (£-77) = a + s T Au a (£). (6.19b) 
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Eliminating the term s T Am (£), we can derive the relations 

(p- q + u 00 (£_- A) -u o(£- m))v o (0 = (p-a)vo(C-A) 

-fe-a)Va(^-M) , (6.20a) 

0-r + w o(£-A)-uoo(£-»7))v a (0 = (p-a)v a (£-A) 

-(r-o)V a (e-»j), (6.20b) 

(r-g + u o(C-?y)- : uoo(C-M))va(C) = (r - a) v a (£ - 7?) 

-(5-o)v„K-A*)- (6.20c) 

We now set p = a then (|6.20a|) and (|6.20b|) become 

p-q + u 00 (Z-\)-uoott-») = -(q-p f P , (6.21a) 



p-r + 5oo(£-A)-u o(£-^) = -(r-p)^_5), (6.21b) 



MO 
P (£-r. 
v P (0 



The combination of (I6.21a[) and (|6.21b|) gives 

p - q + M o(g — A) — u 00 (g - /x) _ p - q + M 00 (£ - A - rj) -Uqo(£ - fi - rf) 
p-r + w o(£ - A) - uoo(C - 17) p - r + 2 00 (£ - A - (i) - %o{£, -»?-/*) 

which is the "lattice KP equation" , [11] . cf. also |10] , 



(6.22) 



From the definition of the function v p (£) in (|6.3j) . (]6.21aj) and (|6.21b[) can be written in 
terms of the r-function 

P-g + «oo(e-A)-u o(e-M) = -(g-p) r( lT A T ) ~^ (0 u > ( 6 - 23a ) 

p-r + u 00 (e-A)-Hoo(C-ry) = -(r - p) r( !~ A ~ v) ^ . (6.23b) 
From (|6.20cp , if we set r = a we also have 

r - q + «oo(£ - A) - u 00 (e - M) = -(9 - r f^l'/'^*® . (6.24) 

US - W ^(4 - ??) 

The combination of (I6.23a|) (|6.23b|) (EOIl) yields 

(p - q)r(£ - A - /i)r(£ - 77) + (r - p)t(£ - A - r?)f (£ - n) 

+(r - g)% - /x - »»)r(C - A) = , (6.25) 

which is actually the bilinear lattice KP equation, (originally coined DAGTE, cf. [5]). 

By introducing the variable r in (|6.ip , we derive the linear and nonlinear relations leading 
to a family of the discrete KP equations (see also [T]). In contrast to the connection estab- 
lished here between the discrete-time RS system and the fully discrete lattice KP equations, 
the rational CM system treated in [21] is connected to a semi-discrete KP equation which 
can be obtained by performing a continuum limit on one of the discrete variables. The con- 
nection between the RS system and solitons has also been discussed in [18]. In the Section, 
we present only the discrete versions of the KP equations. For the continuum limit of these 
equations, it has been studied in [T3]. In the next Section we will how this continuum limit 
corresponds to the non-relativistic limit of the RS system, which after all is viewed as the 
relativistic variant of the CM system. 

7 Non-relativistic limit 

In order to perform the non-relativistic limit A — > 0, it turns out that the Ansatz for the 
Lax pair in Section [2] is a bit too restrictive. For that reason, we will discuss here a slight 
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generalisation of the model by choosing different relativistic parameters for each discrete- 
time flows. The Lax matrices for the " ~ " direction are given by (|2.1bl) and the Lax matrices 
for the " ~ " direction are given by (|2.16bl) with 

N , , N h b 

where /j is the relativistic parameter corresponding to " " direction. Furthermore, we now 
introduce the third discrete-time flow " ~ " which the Lax matrices are given by 

F = V — /'.',,, , and # = V —M—Ey , (7.2) 

f-r 1 x l —Xj+rj ^ Xi - Xj + 77 

where 77 is the relativistic parameter corresponding to " ~ " direction and /, can be deter- 
mined by using the same method of what we did for hi and ki (see Section [5]) . 

Using the same steps as we did in Section [3J we find the system of discrete equations 

P TT (gj ~ x j + A ) _ TT ( x i " x j)( x i — £j + A ) /y o \ 

£ / = i (1* - x j - A ) ~ / = i fat - Sj'X^ - % - A ) ' 

q -i-r (Xj - Xj + fl) _ -i-r (jj - Xj ) (Xj - + /x) . ^ 

5 7=1 (- Ti ~ ^ ~ jJi ( Xi ~ - % - m) ' 

r -pr (gj - + 77) = yX- (,T t - - Xj +7?) , ? g , 

r 7=1 (- Ti ~ ^ ~ ^ 7=1 ~ -^j-fi)' 

where r is the lattice parameter for the " ~ " direction and (|7.3bl) . (|7.3c|) and (|7.3c[) are the 
equations of the motion corresponding to tilde, hat and bar directions with their relativistic 
parameters respectively. The matrix Y of the exact solution fl2.25a|) now takes the form 

Y(n, m, h) = {pi + A)- n (qI + A)~ m (rJ + A)- l Y(0, 0)(pl + A) n {ql + A) m (rl + A) 1 
-np\{pl + A)^ 1 - mqn(ql + A) -1 - lrn(rl + A) -1 , (7.4) 

where I is the discrete-time variable corresponding to the " ~ " direction. 

If we now take the limit on the relativistic parameters such that 71 — > 00, p — > 00 and 
A^O and we define An = — £ and n/p —> t. Then (|7.4[) becomes 

Yfat,m,h) = e- At {qI + A)- m (rI + A)- l [Y(0,Q)+£I}(qI + A) m {rI + A) l e At 
—mqfi(ql + A) -1 — lrr)(rl + A) -1 , 
^Y(£,m,h) = {qI + A)- m (rI + A)- l [Y(0,0)+£I}(qI + A) m (rI + A) 1 

—mqn(ql + A) -1 - lrrj(rl + A) -1 . (7.5) 

(|7.5[) results from the fact that the exponential can move through annihilating each other 
and it is a function of one continuous £ and two discrete-time m, I variables. In fact, the 
limit A — ► coincides with what is some time called 11 the straight continuum limit" of the 
corresponding lattice equation [TT] . 

If we now choose /i = 1/q and 77 = 1/r (I7.5P becomes 

Y(£,m,h) = (ql + A)~ m (rl + A)~ l [Y(Q, 0) + £I](qI + A) m (rl + A) 1 — m(ql + A) -1 
-/(rJ + A)" 1 , (7.6) 
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which is exactly the exact solution for discrete-time Calogero-Moser system [24] [13] . 



From the definition of (|6 . 1 7[) together with (|A.5|i . we can express 



N 



£ - Xi(n,m) 



(7.7) 



and from the definition of (16.31) . we can also express 



<i> = « o (0 = i - £ 



£ - Xi(n, m) 



(7.8) 



where b = /i T (a + Lo)^ 1 h. Eqs (|7.7|) and (|7.8I) form the solution for the lattice KP equation 
(I6.22j) . In the non-relativistic limit we can show that h\ — > 1, cf. [TS], and hence in this limit 
Eq. (I7.7[) becomes 



which is the pole-solution for the semi-discrete KP equation as considered in [T3J [21] • 

The conclusion we can draw from this analysis is that the non-relativistic limit performed 
here provides an explanation as to why the discete-time Calogero-Moser system is derived 
from the pole reduction of the semi-discrete KP equation, while the Ruijsenaars-Schncidcr 
system connects to the fully discrete KP equation, as we have shown in Section [6j In fact, 
the non-relativistic limit engenders in a natural way a continuum limit in at least one of the 
discrete variables of the system. 

8 Discussion 

In this paper we have studied the Lagrangian structure for the Ruijsenaars- Schneider system, 
and shown that similarly to the Calogero-Moser system, which was treated in [21] , it possesses 
a Lagrangian 1-form structure, both on the discrete-time level as well as in the continuous- 
time case. Thus, this is the second example of a system of ODEs which exhibits a Lagrangian 
multi-form structure in the sense of [7] but in a lower-dimensional situation. The present 
example is important, because in contrast to the CM case where the Lagrane structure is 
closely related to the Lax representation (and hence inherit the closure relation from the 
zero-curvature condition), here the relation between the Lax matrices and the Lagrangians 
are less clear, and the validity of the closure relation is more surprising. Thus, we believe 
that these results seem to confirm once again that these Lagrangian form structures are 
fundamental and ubiquitous among integrable systems. 

It is well known that the classical RS system is Liouville-integrable in the continuous-time 
case, [213 121] and formally so in the discrete-time case [TS] HB] as well. However, a note of 
warning is in order here, pointing to the fact that the discrete-time RS system is fundamen- 
tally a complex-valued system of equations, where the solutions may have trajectories which 
even if the initial values are chosen to be real may wander off in the complex plane in finite 
time. In such cases, related to the non-selfadjointness of the underlying eigenvalue problems 
of the Lax pair, it is difficult to make any rigorous statements about the global behavior of 
the solutions. Nevertheless there are perfectly good reasons to study such systems with non- 
physical parameters in a complexified phase space, in particular in connection with problems 
arising from algebraic geometry, cf. e.g. |22j . In the present paper we are interested in these 
models on a more formal level, where the aim is to establish the emergence of a novel struc- 
ture (namely that of Lagrangian 1-forms) and to use these models primarily as toy models 
in order to study how the corresponding variational principles should be implemented. 

With regard to the continuous-time model, the Lagrangian 1-form structure in that case 
had to be established in a rather indirect way, namely by performing systematic limits on 
Lagrangians of the discrete-time system. We have already pointed out that establishing these 




(7.9) 
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Lagrange structures by Legendre transformation from the known Hamiltonians of the model 
is not possible, because it is not a priori known how these Hamiltonian flows are embedded in 
a coherent structure, such that we get acceptable Lagrangian components of the 1-form. In 
this sense the discrete-time model can be viewed as a generating object for such Lagrangians 
for the continuous-time model. Once those Lagrangians have been found, in fact for a model 
with unphysical parameters, we expect that it would be perfectly possible to back-engineer 
the results to find the proper Lagrangian structure for the continuous model in the physical 
regime of the parameters. We are currently investigating this possibility. 

In our view, the importance of this new Lagrangian form structure is that, although 
a formal proof of this assertion has to be found, it seems clear that this structure is the 
manifestation of the multidimensional consistency of the equations in the sense of of the 
papers [12l [2] . It answers the problem of how to find a single Lagrangian framework for a 
situation where we have a multitude of compatible equations imposed on one and the same 
(possibly vector- valued) function of many independent variables. In the case of ODEs, which 
is the case we are dealing with in the present paper, the structure is that of a Lagrangian 
1-form describing systems commuting flows in many time- variables (as many as the number 
of degrees of freedom of the system). The Lagrangians, arising as components of the 1-form, 
should obviously have a very specific forms, in order for the closure relation to hold subject 
to the equations of the motion. In fact, such admissable Lagrangians are themselves in a 
sense solutions of the variational principle underlying the structure, namely the one where 
we consider variations not only with regard to the dependent variables, but with regard to 
the independent variables as well - in other words where we vary the underlying geometries 
in the space of independent variables. This forms in our view poses a new paradigm in 
variational calculus, [3], and possibly a new principle of fundamental physics. That theories 
which exhibit such structures always correspond to integrable systems (i.e. integrable in 
other commonly used senses of the term, e.g. through inverse scattering, existence of Lax 
pairs etc.) is challenging question that remains to be answered. 



A The construction of the exact solution 

In this Appendix we review the construction of the exact solution for the RS system. The 
basic relations following from the Lax pair (|2. lb|) together with the definitions (|2.2j) . lead to 

XM + XM — M Q X = lih T , (A.l) 
XM Q + XL - L a X = hh T , (A.2) 

where we introduce X = X)t=i x i^a- O n t ne other hand, from the Lax equation (|2.2[) . we 
obtain the relations 

L M = M L , (A.3a) 
L h - M h = -ph , (A.3b) 
h T L - h T M = -ph T . (A.3c) 

We now introduce another transformation 

L = UoAUq 1 , and M = U a U^ , (A.4) 

where Uq is an invertible N x N matrix, and where the matrix A is constant: A = A, as a 
consequence of (|A.3b[) . Then introducing 

Y = Uq 1 XUq , r = Uq 1 • h , s T = h-U , (A.5) 

we obtain from (|A.3j) and (|A.4j) . 

(pi + A) • r = r , s T ■ (pi + A) = s T , (A.6) 
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where I is the unit matrix, as well as from (IA.1I) and (IA.2I) . we have 

X + Y — Y = rs T , (A.7) 
AA + [Y, A] = rs T . (A.8) 

Eliminating the dyadic rs T from ()A.7|) by making use of (|A.6[) . we find the linear equation 

Y = ( p I + A)- 1 Y(pI + A)--^- , (A.9) 

pi + A 

which can be immediately solved to give 

Y(n, m) = (pi + A)"" ( Y(0, m) - ™ pX A ) (pi + A)™ , (A.10) 
V pI + AJ 

subject to the constraint on the initial value matrix 

[Y(0,m),A] = AA + rank 1 . (A.ll) 

A similar analysis can be applied to create the solution associated with the . 

Conversely, we can start from a given N x N diagonal matrix A with distinct entries, 
and an initial value matrix 1^(0,0) subject to the condition that 

[Y(0, 0), A] = AA + rankl , (A.12) 

where [ , ] represents the matrix commutator bracket. Let J7 _1 (0, 0) be the matrix that 
diagonalized Y(0,0), i.e., such that 

F(0,0) = J7 _1 (0,0)X(0,0)J7(0,0)) , X(0, 0) = diag(xi(0, 0), . . . ,x N (0, 0)) . (A.13) 

If the eigenvalues of 1^(0, 0) are distinct (which we can take as an assumption on the initial 
condition) then (7 _1 (0,0) is determined up to multiplication from the right by a diagonal 
matrix times a permutation matrix of the columns. (Fixing an ordering of the eigenvalues 
Xi(0, 0), C/ _1 (0,0) unique only up to multiplication by a diagonal matrix from the right). 
We can fix C/ _1 (0, 0) up to an overall multiplicative factor by demanding that 

[Y(Q, 0) , A] = AA - r(0, 0) s T (0, 0) . (A.14) 

Next, we consider the matrix function given by 

Y(n,m) = (pI + A)- n (qI + A)- m Y(0,0)(pI + A) n (qI + A) m 

-np\(pl + A)- 1 - mq\(ql + A)" 1 (A.15) 

Let U(n,m) be the matrix diagonalizing Y(n,m) by an appropriate choice of an overall 
factor (as a function of n and m) this matrix can be fixed such that it obeys: 

r(n,m) = (pi + A)- n (qI + A)-"Y(0, 0) , and s T (n,m) = s T (0,0) (pi + A) (ql + A) , 

(A.16) 

and 

[Y(n, to) , A] = AA — r(n, m) s T (n, to) . (A. 17) 

From the expression (|A. 15|) we can derive the relations 

(pl + A)Y Y (pl + A) = - P X , (A.18a) 
(ql + A)Y-Y (ql + A) = -qX, (A.18b) 

with the usual notation for the shifts in n and to over one unit. Together with the relation 
(|A.17|) this subsequently yields: 

Y-Y = -pX + rs T , Y - Y = -qX + rs T . (A.19) 
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Reversing these relations by rewriting them in terms of 

X{n,m) = U{n,m)Y(n,m)U- 1 {n,m) (A.20) 
and now defining the Lax matrices by 

L:=UAU l , K-UAU 1 , (A.21) 



together with 
we recover the relations: 



and 



M;=UU~ l , N ;= UU^ 1 , (A.22) 

[X , L] = XL + hh T , (A.23a) 
[X , K] = XK + hh T , (A.23b) 

XM-MX = -XM + hh T , (A.24a) 
XN-NX = -XN + hh T , (A.24b) 

which determine the matrices M and N as functions of the Xi(n,m) as well as the off- 
diagonal parts of the matrices L and K. 
Furthermore, from (|A.18[) we obtain 

LXM - MXL = -phh T , (A.25a) 
KXN - NXK = -qhh T , (A.25b) 
which, when combined with the relations of (|A.23[) . yield 

(LM - ML) (X + X) + (Zh - Mh \ h T = -phh T , (A.26a) 

(KN - NKj (X + X) + (Wi - Nh\ h T = -qhh T . (A.26b) 

On the other hand, using the relations (IA.24I) we also obtain 

(X + X\ (LM - ML) +h(hL- h T M) = -phh T , (A.27a) 

(x + A) (KN - NKj + h (hK - h T N) = -qhh T . (A.27b) 

From the relations (|A.26[) and (|A.27p it follows that the Lax equations hold and their form 
is determined up to the diagonal part of the matrices L and K . 

B The compatibility between the Lax matrices 

In this Appendix we will show the explicit computation of the compatibility between Lax 
matrices. 

B.l The compatibility between L K and K K 

Consider the relation 

L K K K = K K L K , (B.l) 

yielding 

N f 1 



K K + 2X + Xi — xi Xi — Xj + A 



= fc 2 y ( -i + — - — . (b.2) 

^— ' J V k k + 2A + Xi — xi Xi — x ; + X J 

3=1 
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(|B.2jl leads to the condition 

N N 
3=1 3=1 

Introducing j3 = Y^j=i hj/12j=i j we nave 

ft? - ^ ■ 

B.2 The compatibility between L K and iV K 

consider the relation 



L K N K = N K L K , 



yielding 

N 



k k + 2A + Xi — xi xt — Xj + X xi — Xj + X 

3=1 v j j 

" '• i ii 



n + 2X + Xi — xi Xi — Xj + X Xi — xa + 

j=l \ 

The terms with the k lead to the condition 

hi = ak? , 

where a = ^j/ Sjli we have 

ft? = afc? • 

Comparing (IB.8|) with (|B.4[) . we have a = [} which implies Ylj=i hj = ^ 
The remaining terms in (|B.6|) lead to the identity 

x / afc? afc? \ JL / afc? afc 2 N 



y ^—^i —^i — =y 



J J 



. Xi — Xi ■■ + X Xi — Xj + X t— f \ Xj — arr + A x< — xi + A 

j=l \ J J / 3 = 1 \ j 

Using the Lagrange interpolation formula, we have 

p _ Uf=i fa - x 3 - X)(xj - Xj + A) 

rij^feiC 3 '* _ ^j) rij=i fa _ ^i) 
2 _ n^Li fa - %• - x )( x i - x 3 + A ) 

from which we obtain the equations of motion (|2.19l) . 

B.3 The compatibility between K K and M K 

consider the relation 



K K M K = M K K K , 



yielding 

JV 



m £ Vj ( i — _u — + + ^-^) 

J \ K K + ZA + Xi — Xl Xi — Xj + X Xi — Xj + A J 

~ - f 1 



1 1 1 



K fx ~\~ 2 A ~\~ Xi X\ Xi ~~~Xj -f~ A Xi Xj ~\~ 
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The terms with the k lead to the condition 



where /i = Y,j=i tf/ YljLi &f > we nave 



Comparing (IB.15|) with (|B.4|) . we have M = /3 which implies Ylj=i &f = 2~Zj= 
The remaining terms in (|B.6I) lead to the identity 



1,2 
i K i- 



y 



E 



- a;; + A - + A 



Using the Lagrange interpolation formula, we have 



m, 2 



m, 2 



- X)(xi — Xj + A) 



from which we obtain the equations of motion (|2.12l) . 



C Examples 



(B.14) 



(B.15) 



(B.16) 



(B.17) 
(B.18) 



m 




(n ,m ) 



N«- 



r' 



mi 



N - 1 = Ni 

(b) 



(No, w ) 



Figure 2: The effect of changing variables on the discrete curve I. 

In this Appendix, we will investigate how to derive the discrete Euler-Lagrange equation 
from the variational principle. For general discrete curves it is cumbersome to implement the 
variational principle because of the notation it would require. We will, however, demonstrate 
how the principle works for a few simple cases: (a) the curve shown in Fig. (b) the 
curve shown in Fig. ([3]). 

Case(a): The curve shown in Fig. (0): We now introduce a new variable N = n + m 
(which will play an important role in the next section) together with the change of notation 



x(n,m) i-> cc(N,m), x := £c(N + 1, m) and x := x(N, m + 1) 
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and so we work with the curve given in Fig. (JSt)- The action evaluated on this curve can 
be written in the form 

mi — l mi — 1 

S[x;T'] = -^(M)(x(N a -l,m),x(N Q ,m))+ ^ jSf {m) (x(N - l,m),x(N 0> m+l))(C.l) 



where 



N N 



N N 



'•.;=i 



(C.2) 



N N 



i,j=l 



. N N 

^ X ~ y J + A ) - ln qVP\ ^2( x i ~ Vi) > 



(C.3) 



The minus sign in f|Cl[) indicates the reverse direction of the Lagrangian L(n) along the 
horizontal links. Performing the variation x t— > x + Sx, we have 



SS = Q = 

■m i — l 

X 

m—m,Q 



&S? (N) (a;(No - l,m),x(H ,m)) 



Sx(No,m) 



8x(Nq, to) 

ajSf (N) (aj(No- l,m),x(N ,m)) 



8x(Nq — 1, m) 
mi_1 /a=Sf (m) (a;(No - 1, m), s(N , m + 1) 



&c(N - l,m) 



+ ? V 



aa;(N , m + 1) 

9jSf (m) (a(N - l,m),a;(No,m + l) 
+ 0x(N o -l,m) 



6x(N - l.m) 



(C.4) 



We now obtain the Euler-Lagrange equations 

9if( N) (a;(No - l,m),x(N ,m)) 9JSf( m )(sc(N - l,m- l),x(N 0) m) 



+ 



<9:c(IMo,m) <9:c(IMo,to) 
aj2f {N )(x(N -l,m),x(No,m)) flJ2f (m )(x(N - l,m),x(N ,m + l) 



which produce 
In 



<9a:(No - 1, to) dx(N - 1, m) 

aJi(No,m) - aCj(N - l,m) 



, (C.5a) 
= , (C.5b) 



In 



_P_ 



N 



N 
3 = 1 



Xi(No, to) — Xj(No — l,m) + A 

Xi(N ,m) - xj(No - 1,™ - 1) + A 



£Ci(N ,m) - Xj(N - l,m- 1) 
Xi(N Q -l,m) - Xj(N ,m+ 1) 



a;, (N — l,m) — Xj(Uo, m + 1) — A 

Xi(N — 1, to) — x 3 -(No, to) — A 



Xi(H - l,m) - Xj-(N ,m) 
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which are equivalent to (I2.24a|) and (|2.24b[) . respectively. 




(n ,m ) 



mi 



(N' OJ mo) N' 1 = N(, + 1 
(b) 



N' 



Figure 3: The effect of changing variables on the discrete curve II. 

Case(b): The curve shown in Fig. ([3]): Introducing the variable N' = n — m, the corre- 
sponding curve is given in Fig. ([3b). The action evaluated on the curve T' reads 



S[x;T'} - ^(wMK,m),x(N' + l,m)) 

m—rriQ 

+ £' (m) {x(N' + l,m),x(N' ,m + l)), 



(C.7) 



where 



N N 

3? m (x,y) = {f(xi-y 3 ) - fixt-yj ~ \)) - - f( yi -y.j + X) 



N N 



1=1 



N N 
■&(m) (x, y) = (f(Xi -Vj)~ f(Xi - Vj - A)) - - Y f( x i ~ x j + A ) 



N N 

, X! ~ y j + A ) - ln i\fp\ Y( Xi ~ y$ ' 



i=l 



(C.8) 



(C.9) 
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Performing the variation x i— > x + 5a;, we have 



(55 = = 

"" 1 /SJSf (N0 (x(N{,,m),aj(N{, + l,m) 



E 



m-mfl 



dx(l%,m) 



5a;(Ng,m) 



djy (N ,)(x(N(,,m),x(N{, + l,m) 
dx{N' + l,m) 



fe(Np + l,m) 



mi —1 

E 



aSf (m) (a;(N^ + 1, m), x(N(„ m + 1] 



5x(N' ,m+ 1) 



dx(]\l£,,m + l) 

dJtf {m) (x(N' + 1, m), a(N(,, m+1) 
+ 0x(l% + l,m) 



*»(Nj, + l,m; 



We now obtain the Euler-Lagrange equations 

dJ2? (N0 (x(N' , m), x(l% + 1, m) aj£f (ro) (x(N' + 1, m - 1), x(N' , m) 



+ 



dx(N(,,m) 9x(l%,m + l) 

d^ NI) (x(N' ,m),x(N' + l,m) 0jS? (m) (x(N& + l,m),x(l%,m + 1) 



dx(N(, + l,m) 



which produce 
In 



?/3 



In 



iL 

?/3 



E(* 

3=1 



3=1 



<9x(l% + l,m) 



Xi(N' ,m) - Xj(H' - l,m) 



= 0, 
= 



(CIO) 

(C.lla) 
(C.llb) 



Xi(N(),m) - Xj(N() - l,m) + X 

Xi(%,m) - Xj(Nq - 1, to - 1) + A 



In 



x i (N' ,m)-x j (N' -l,m-l) 
Xi(K a - l,m) - x J (N / ,m+ 1) 



Xi(N^ - l,m) - x^N^m + l) - A 

sc,(Nq — l,m) — X i j(Ng, to) — A 



In 



Xi(N^ - l,m) - Xj(N(,,m) 



which are equivalent to (|2.24al) and (|2.24bl) . respectively. 

By working with the specific type of curves given in Fig. ^ and Fig. ©, we can perform 
the variational principle with either the new variables (N,m) or (N',to). Even though we 
have these two possibilities, the natural choice is to work with (N,m). This is because the 
variable N' leads to difficulties in performing the skew limit in the next section. 
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